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, Abstract. In this article we review the current status of the stellar population assignment for 

i novae. Observations in the MiUcy Way and in external galaxies point out the existence of two nova 

' populations: fast and bright novae, mainly originated from massive white dwarfs and associated 

with the thin disk/spiral arm stellar population, and slow and faint novae, originated from lighter 
■ white dwarfs and associated with tliick-disk/bulge population. 

(N 

INTRODUCTION 



\ Baade [1],[2] introduced the concept that different kinds of stellar populations have 

different spatial distribution within galaxies (see also Oort [3]). We can take advantage of 
^ . this notion to find out useful hints about the population assignments of the progenitors of 

O I novae. Due to their luminosity, My ^ — 9 novae are particularly suited for this purpose 

fSj ' because they can be easily identified both in the Milky Way and in external galaxies. 

o ■ 

NOVA POPULATION IN THE MILKY WAY 

O 

^ • Historical data on galactic novae (e.g. McLaughlin [4], [5], [6] and Payne- Gaposchkin 

■ [7]) have received discrepant interpretations. Kukarkin [8], Kopylov [9] and Plaut [10] 
pointed out the existence of a concentration of novae towards the galactic plane and the 

k> ^ galactic center and classified them as belonging to the 'disk population' . Minkowski[l 1], 

■ [12] and [7] showed that the galactic longitudes of novae and planetary nebulae (PNe) 
' have similar distributions and therefore novae, like PNe, belong to Pop II stellar popu- 
lation. Baade [13] assigned novae to Pop II stellar population because of the occurrence 
of a few ones (e.g. T Sco 1860) in very old stellar population systems, such as the Glob- 
ular Clusters. Iwanowska and Burnicki [14] suggested that novae are a mixture of Pop 
I and Pop II objects, and Patterson [15] proposed that novae belong to an 'old disk' 
population. Tomaney and Shafter [16] found that novae belonging to the bulge of M31 
are spectroscopically different from novae observed in the neighborhood of the Sun and 
deduced that galactic novae are mainly 'disk' objects. Different conclusions were drawn 
by Delia Valle and Duerbeck [17] who compared the cumulative distributions of the 
rates of decline for M31, LMC and Milky Way nova populations (see Fig. 1) and found 
that galactic and M31 distributions are indistinguishable, whereas M31 and LMC dis- 
tributions are different at ^ 99% significance level. Since the speed class of a nova 
depends on the mass of the underlying white dwarf (e.g. [18]), systematic differences 
in the distributions of the rates of decline indicate the existence of physical differences 
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FIGURE 1. Cumulative distiibutions of the rates of decline for M31, LMC and Galaxy (from[17]) 



between LMC and M31 nova populations. As most novae in M31 are produced in the 

bulge (CiarduUo et al. [19], Capaccioli et al.[20], Shafter and Irby[21]), one would argue 
that novae in the Milky Way are also mostly bulge objects. 



Disk and Bulge Novae 

The quantitative characterization of the concept of nova populations into two classes 
of objects, i.e. fast and bright 'disk novae' and slow and faint 'bulge novae', has been 
elaborated in the early 90's by Duerbeck [22] and Delia Valle et al.[23],[24], [25],[26]. 
The former demonstrated that nova counts in the Milky Way do not follow an unique 
distribution (Fig. 2), the latter authors showed that the rate of decline (which traces the 
mass of the WD associated with the nova system) correlates with the spatial distributions 
of the novae inside the Milky Way (Fig. 3 and 4). Fig. 2 shows that nova counts 
follow two different trends. Dashed and dotted lines are the predictions from simple 
disk, p(z) = a X Poexp( — |z|/zo), and bulge (p ~ 10^-^) nova population models [with 
Po = 125pc Ay = lmag/kpc"\ My(max)= -9, po = 10"^*^ pc~^ yr"^ and a ^0.1, 
see [15], [17], [27] [28]]. 
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FIGURE 2. Theoretical and observed number of novae. Filled circles: observed rate; dashed line: 
calculated number counts for a disk population; dotted line: expected slope for a bulge population. 
Adapted from [22]. 



A simple consequence of two population distributions is that the rate of decline is 
expected to correlate with the galactic longitude. Fig. 3 compares the distributions of 
t3 for galactic novae in the direction of the galactic anti-center (+90° < / < 270°, 
dotted region) and center ( 90° < / < +90°). The former distribution mainly formed 
by 'disk' novae, peaks at t3 ^ ISJ^ , while the latter one peaks at larger t3 because of 

the contribution of sJow 'bulge' novae which are viewed in the direction of the galactic 
center. 



The relationship between the rate of decline and height above the 
galactic plane for classical novae 

If galactic novae originate from both the bulge and thin disk of the Milky Way, then 
the masses of their WD progenitors are expected to have heights above the galactic 
plane systematically different because of the initial-mass/final-mass relationship for 
WDs (e.g. [29]). This fact can be verified in the following way. Theoretical calculations 
([30], [31], [32]) have established that the strength of the nova outburst is a strong 
function of the mass of the underlying WD. On the other hand, the luminosity of the 
nova at maximum (which is representative of the strength of the outburst) correlates 
with the rate of decline ([33], [34], [35]). Thus, the distribution of the rates of decline 
traces the distribution of the masses of the WDs associated with the nova systems. We 
note that the distance moduli introduced into z = iQl^-'^'im-M+s-A)] ^ ^ (where b is 
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FIGURE 3. The distribution of t3 (the time that the nova takes to fall 3 magnitudes below maximum 
light) in the direction of the galactic anticenter (+90° < I < 270°, shaded region) and center (—90° < I < 
+90°, black). Adapted from [23] 



the galactic latitude) are derived via expansion parallaxes rather than with the maximum 
magnitude vs. rate of decline (=MMRD) relationship to avoid circular arguments. For 
this reason Delia Valle et al. [23] restricted their analysis to a fiducial sample of only 19 
objects coming from Cohen & Rosenthal [36] and Cohen [37] out of about a hundred 
novae for which the rates of decline have been measured with reasonable accuracy. In 
this presentation we were able to augment the original sample by 1 1 new objects from 
[38], [39] and [40]. Fig. 4 shows the existence (at a confidence level of ^ 3o) of a 
relationship between the rate of decline and the height of the novae above the galactic 
plane, log z = — 0.58(±0.18) x (log 100 x v^) +2.7(±0.2), which can be expressed as 
log z = -1.9(±0.5) X M^^jj/Mq + 3.7(±0.6), after recalling that oc log (from 
[41]) and Mg log t2 (from [35]). The trends illustrated in Fig. 4 and 5 indicate that the 
fastest nova systems, which contain the most massive WDs, are concentrated close to the 
galactic disk. It is difficult to understand how such distributions could be due to selection 
effects, since there is no obvious mechanism to prevent the discovery of fast/bright novae 
at high z (although we cannot exclude that some slow/faint nova at small z can be heavily 
absorbed and overlooked). Fig. 6 reports the frequency distribution of the rates of decline 
(in terms of log t2) for the fiducial sample of novae. The distribution is bimodal and 
shows that 'disk' novae, i.e. objects characterized hy z ^ 150pc (shaded region), are 
mostly 'fast' whereas 'bulge' novae are mostly 'slow'. The former distribution peaks at 
My = —^-T^q'I and the latter at = — 7.2^^ The existence of such a bimodality was 
pointed out by Arp [42] also for M31 novae (see his Fig. 36), but this fact was neglected 
afterwards. 
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FIGURE 4. The relationship between height above the galactic plane vs. rate of decline, ~ 2/t2. t2 
is the time that the nova takes to decrease its brightness by 2 magnitudes from maximum. The sizes of the 
errorbars in log z axis are comparable to the sizes of the dots. 




FIGURE 5. The relationship between height above the galactic plane vs. mass of the underlying WD 
(in solar masses). The sizes of the errorbars in log z axis are comparable to the sizes of the dots. 



The spectroscopic differences between disk and thick-disk/bulge 

Novae 



Williams [43] after studying about two dozen of galactic novae (and a few objects 
belonging to the LMC) concluded that novae can be broadly divided into two spectro- 
scopic classes of objects: the Fe II and He/N novae. The former are characterized by 
slow spectroscopic evolution with expansion velocities ^ 2500 km/s (FWZI) and the 
Fe n lines as the strongest non-Balmer lines in the early emission spectrum. The lat- 
ter are fast spectroscopically evolving novae, characterized by high expansion velocity 
ejecta ^ 2500 km/s (FWZI) with He and N lines being the strongest non Balmer lines 
in the emission spectrum near maximum. Hybrid objects (e.g. V1500 Cyg) that evolve 
from Fe II to He/N are classified as Fell-b (b=broad) and are physically related to the 
He/N rather than to Fell class. In Fig. 7 we have plotted the frequency distribution of the 
heights above the galactic plane of the novae of the fiducial sample after being classified 
according the Williams' criteria. The top and bottom panels give the distributions of no- 
vae which have been classified as He/N (-i-Fell-b) and Fe II. The histograms show that 
novae belonging to the He/N class tend to concentrate close to the Galactic plane with a 
typical scale height ^150 pc, whereas Fell novae are distributed more homogeneously 
up to z ~ 1000 pc and beyond. A K-S test on the data shows that the two distributions are 
different at ^ 95% level. Fig. 7 indicates that the objects previously classified as 'disk' 
and 'bulge/thick-disk' novae tend to correspond to the spectroscopic classes introduced 
by Williams, indeed about 70% of fast and bright novae belong to the He/N (or Fell-b) 
class, while the slow and faint ones form the main bulk of Fe II class. To explain this be- 
havior one should consider the following. The more massive the WD (for a given M and 
T^y^) the smaller is the mass of the accreted envelope (in view of Miacc °^ ^wd/^wd^ 
the more violent is the outburst (i.e. shorter t3 and higher expansion velocities) and the 
larger the fraction ^shock/ ^^wind where Am^^^^^ is the shell mass ejected at the maxi- 
mum and Am^^.-^^^ is the fraction of the shell mass ejected in the subsequent continuous 
optically-thick wind phase. Since He/N spectra are formed in the shell ejected at the out- 
burst maximum (as one can infer from the top-flatted profiles of the emission lines, for 
example) it is very likely that He/N novae are generally associated with massive WDs. If 
this interpretation is correct the distributions reported in Fig. 7 are simple consequences 
of the trends reported in Fig. 4 and 5. 



NOVA POPULATIONS IN EXTERNAL GALAXIES 

The use of galactic data to establish the stellar population where novae originate has 
been often questioned because of observational bias, mainly due to both interstellar 
absorption in the galactic disk and our position within the Galaxy. These effects can be 
largely minimized by studying the nova populations in external galaxies, particularly: 
1) their spatial distribution; b) the maximum magnitude vs. rate of decline relationship; 
c) differences (if any) in the nova rates exhibited by galaxies of different Hubble type. 
The available data are summarized in Table. 1. Col. 1 gives the galaxy identification; 
col. 2 the nova rate; col. 3 the total B mag of the galaxy corrected for background and 







1 2 
log 12 



3 



FIGURE 6. Frequency distribution of the rates of decline for the novae of the fiducial sample. The 
shaded region represents novae with z ^ 150pc. 



TABLE 1. Absolute and normalized nova rates 



Galaxy 


Novae/yr 


^tot 


(B-K) 


(m-M) 




T 


LMC 


2.5 ±0.5 


0.57 


2.74 


18.58±0.1 


5.1±1 


9 


SMC 


0.7 ±0.2 


2.28 


2.71 


19.00±0.1 


4.8 ±1.5 


9 


M33 


4.6 ±0.9 


5.75 


2.87 


24.64 ±0.2 


3.7 ±0.9 


6 


MlOl 


12±4 


8.26 


3.24 


29.35 ±0.2 


0.9±0.3 


6 


M51 


18±7 


8.41 


3.43 


29.60 ±0.2 


1±0.4 


4 


M31 


29±4 


3.51 


3.85 


24.42 ±0.2 


1.5 ±0.2 


3 


M81 


24 ±8 


7.39 


3.99 


27.80 ±0.2 


1.7 ±0.6 


2 


N5128 


28±7 


6.32 


3.38 


27.80 ±0.2 


3.3±0.8 


-2 


N1316 


130±40 


9.20 


4.15 


31.50±0.3 


1.7 ±0.5 


-2 


M87 


91 ±34 


9.49 


4.17 


30.90 ±0.2 


2.1 ±0.8 


-A 


VirgoEs 


160 ±57 


9.46 


4.26 


31.35±0.2 


2.2±0.8 


-A 



foreground absorption; col. 4 the color (B-K); col. 5 the adopted distance modulus; col. 
6 the normalized nova rate, i.e. the nova rate per year/1 O^'^L^q; col. 7 the Hubble type. 
Data come from [24] and [44]. Note that to be consistent with the results of MACHO[45] 
and EROS [46] teams, who found 4 novae in the SMC after one year of monitoring to 
search for microlensing events, the nova rate in this system has been estimated, on the 
basis of simple statistical arguments, to be ~ twice as large as quoted in [44]. 
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FIGURE 7. Frequency distribution of the height above the galactic plane for He/N (and Fell-b) novae 
(top panel) and Fell novae (bottom panel). Adapted from [26]. 



Due to its closeness, the M31 nova population has been the only one extensively 
studied. However a simple skimming through past literature reveals the existence of 
different ideas on its nova population assignment: disk populations (Arp [42], Rosino 
[47]);halo population (Wenzel and Meinunger [48]); mainly bulge population (CiarduUo 
et al. [19], Capaccioli et al. [20] see their Fig. 7); mostly disk population (Hatano et al. 



It is very likely that a fraction of M3 1 novae considered by [19] and [20] as belonging 
to the bulge are in fact (according to [49] see their Fig. 1 and 4) physically related to 
the disk and their allotment to the bulge was a simple consequence of neglecting the 
geometrical projection effect. However Shafter & Irby [21] have reappraised the Hatano 
et al's conclusions, after comparing the spatial distribution of novae with the background 
M3 1 light (see their Fig. 4). They confirm most M3 1 novae to be originated in the bulge 
(about 70% up to a minimum of 50%). Concerning this last point a word of caution 
seems in order as long as the effects of extinction on nova detections in the M3 1 disk 
will be fully quantified. 

M33 is an almost bulgeless galaxy [50] and therefore its nova production necessarily 
originates in the disk. Particularly Fig. 2 of [24] shows that most novae appear superim- 
posed on the arms of the parent galaxy. Similar arguments hold for the LMC. 

The analysis of 78 plates obtained in the 1950-55 Palomar campaign for the discovery 
of novae in M81 has been recently published by Shara, Sandage and Zurek [51]. These 
authors find that the spatial distribution of novae in M81 is fully consistent with the two 
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FIGURE 8. The MMRD relationship for LMC, M31 and Virgo nova populations. The solid Une 
indicates the best fit (from [35]). Upper and lower dashed curves are located at ±3(7. 



nova populations hypothesis (see their Fig. 1). Particularly they estimate that the fraction 
of novae belonging to the M81 bulge is not larger than ~ 60%. 



The Maximum Magnitude vs. Rate of Decline Relationship 

Fig.8 reports the maximum magnitude vs. rate of decline relationship for LMC, M31 
and Virgo novae. A simple glance indicates that the nova production in the LMC is 
clearly biased towards fast and bright novae, whereas the M3 1 nova population exhibits 
a prominent 'slow' component. By the same token, admittedly on the basis of a scanty 
statistic, the trend exhibited by novae discovered in Virgo (Pritchet & van den Bergh 
[52]). It is not obvious to explain this behavior in terms of an observational bias: indeed 
the brightest novae are detected in the nearest galaxy (LMC) and would be missed in 
the more distant ones (Virgo). On the other hand, the differences in the MMRDs find a 
simple explanation in the framework of the two nova populations scenario: nova systems 
in disk dominated galaxies are associated with more massive WDs, then resulting in 
faster and intrinsically brighter nova events. 
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FIGURE 9. The nova rate per unit of K luminosity as a function of the Hubble type of the parent galaxy. 
The error bars reflect only the poissonian errors affecting the single nova rates. 

Is the nova rate depending on the Hubble type of the parent galaxy? 

Delia Valle et al. [24] suggested that there may exist a systematic difference in the 
nova rate per unit of K (and H) luminosity in galaxies of different Hubble types (see their 
Fig. 3 and 4). These authors found that galaxies of late Hubble types are more prolific 
nova producers than early-type ones by a factor ~ 3. This overproduction of novae could 
be the consequence of "selection effects" on the nova frequency (Truran & Livio [53], 
Ritter et al. [54]). For example, if we use the results of [54], we obtain for the ratio of the 
nova outbursts density, P(^^«'S'^)oMfte«f/P(^"^^^)oM?Z)Mwr ~ 4 in good agreement with the 
observations. A possible interpretation is that nova systems in disk dominated galaxies 
result in more frequent nova events due to the shorter nova recurrence time associated 
with massive WDs (see [55]). However this finding has been put into question by [44] 
(see also [56]) on the basis of nova rates on M51 and MlOl. Fig. 9 shows the trend of the 
normalized nova rate vs. (B-K) as derived from Tab. 1. With the exception of MlOl 
and M5 1 (for which the nova rates reported by [44] may be lower limits) the normalized 
nova rate increases from ^ 2 novae yr^^/lO^'^L^Q for early types to about ~ 5 
for late types. 



Both observations of galactic and extragalactic novae converge toward the idea that 
novae are drawn from two different stellar populations. In particular: 



CONCLUSIONS 



1) analysis of the nova counts and rates of decline inside the Milky Way has allowed 
to derive the notion of disk and bulge/thick disk nova populations. The typical disk 
nova is a fast evolving object whose lightcurve exhibits a bright peak at maximum 
M^, > - 9 (t2 < 13'' or t3 ^ 20'^) a smooth early decline and belongs to the He/N 
(or Fe lib) spectroscopic class. The progenitor is preferentially located at small heights 
above the galactic plane ( ^ 150pc) and since it is related to (a relatively old) Pop 

I stellar population, the associate WD is rather massive, Mjy^ ^ 1 M . . The typical 

bulge/thick disk nova is a slow evolving object whose lightcurve exhibits a fainter peak 
at maximum. My ~ —7.2 (t2 ^ 13^ or t3 ^ 20^), often double maxima, dust formation, 

maximum standstill and belongs to the Fe II spectroscopic class. The progenitors extend 
up to 1000 pc from the galactic plane and are likely to be related to a Pop II stellar 
population of the galactic thick-disk/bulge and are therefore associated (on average) 
with less massive WDs, M^^ <, IMq. 

2) analysis of the MMRD relationship for LMC, M31 and Virgo novae confirm the 
existence of systematic differences in the distributions of the rates of decline of the 
respective nova populations. 80% of novae in the LMC are bright and fast and therefore 
associated (on average) with massive WDs while fast novae in M31 are ^25%. The 
two distributions are significantly different (a K-S test gives ^ 99%). The analysis of 
the spatial distributions suggests that nova populations in M31 and M81 are a mixture 
of disk and bulge novae ( ~ 30% — 70% and 40% — 60% respectively). Novae in the 
LMC and M33 originate from disk population whereas novae in M87 and NGC 1316 
are related 'a fortiori' to bulge population. 

3) the previous results imply the existence of differences in the age of the progenitors 
of disk and bulge novae. This has been recently proven by Subramanian & Anupama 
[57]. These authors have studied the star formation history of nearby regions around 
LMC novae and show that most of them occur in the bar or close to it. They conclude 
that the parent stellar population of the fast and slow novae are likely to be in the range 

<, 1 — 3 and >, 3 — 10 Gyr respectively. 

4) some of the differences between disk and bulge novae, as described above, are 
expected on the basis of theoretical arguments (see Kolb [58], Starrfield et al. [59], Kato 
[60]). 

5) Analysis of the distributions of the rates of decline as a function of the absolute 
magnitude at maximum of LMC, M31 and Virgo, reveals the existence of a scant 
group ( ^ 5%) of super-bright novae which deviate systematically from the MMRD 
relationship by more than one magnitude [61]. One possible explanation is that a "super- 
nova" explosion might occur at the end of the life of a CV ([62], see also[63]). The 
recent study on Nova LMC 1991 ([64]) suggests that the metallicity may be the driving 
parameter to account for this deviating behavior. 

Future studies 

Whether or not the nova rate depends on the Hubble type of parent galaxy is still 
an open question. If low nova rates for M51 and MlOl will be confirmed a possible 



explanation for the trend reported in Fig. 9 has been pointed out by Yungelson, Livio 
and Tutukov [65] (see also [66]). These authors noted that small are associated with 
high-mass spirals such as M51 and MlOl, while high values of are typical of low 
mass galaxies such as LMC and M33. If the presently observed rates of disk novae in 
late type galaxies are mainly determined by the current SFRs, one expects to observe 
differences in the normalized rates of high- and low-mass spirals. Indeed for the former 
the SFR in the past was several ( ^ 5) times higher than the present one (Gallagher 
et al. [67]). For the latter the SFR was nearly constant over the galaxy lifetime (e.g. 
Gavazzi & Scodeggio [68]). Therefore high-mass spirals contain a higher fraction of old 
red stars than low-mass ones and consequently also their IR-luminosity is higher (this 
is actually observed, see col. 4 in Tab.l). As a consequence, since the K luminosity is 
proportional to the mass in old stars, one may expect in systems such as M33 or 
LMC to be higher than that in M51 or MlOl. It is apparent that this issue will be solved 
as soon as new nova rates in external galaxies will be made available. For example, 
preliminary results, based on HST observations (e.g. [69]), may suggest high values of 
nova rates also in elliptical galaxies. We note that nova survey in extragalactic systems, 
despite their scientific interest, have not been popular among astronomers (although 
'remarkable' exceptions do exist, e.g. [70]). Probably the main reason for this is the 
unpredictable nature of nova events, which made nova surveys considerably (telescope) 
time consuming. However, the coming into operations of 8- 10m class telescopes should 
change this bias. Delia Valle and Gilmozzi [71] have carried out with VLT a pilot 
programme to discover novae in NGC 1316. They found 4 novae with 3h of observing 
time (about ~ 0.8'^ per nova). A similar programme carried out by [52] with a 4m class 
telescope on Virgo galaxies (which should have comparable nova rates and distance than 
NGC 1316) discovered 9 novae in 56h (~ 6.5'' per nova). This experiment teaches us that 
the 8-lOm class telescopes equipped with larger and more efficient detectors are able to 
improve the payoff per night, in terms of nova detections in galaxies outside the Local 
Group, by on order of magnitude with respect to the previous generation of telescopes. 
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